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Particle delocalization is a common feature of quantum random walks in arbitrary lattices. However, in the 
typical scenario a particle spreads over multiple sites and its evolution is not directly useful for controlled quan¬ 
tum interferometry, as may be required for technological applications. In this paper we devise a strategy to 
perfectly split the wave-packet of an incoming particle into two components, each propagating in opposite di¬ 
rections, which reconstruct the shape of the initial wavefunction after a particular time f. Therefore, a particle in 
a delta-like initial state becomes exactly delocalized between two distant sites after f . We find the mathematical 
conditions to achieve the perfect splitting which are satisfied by viable example Hamiltonians with static site- 
dependent interaction strengths. Our results pave the way for the generation of peculiar many-body interference 
patterns in a many-site atomic chain (like the Hanbury Brown and Twiss and quantum Talbot effects) as well 
as for the distribution of entanglement between remote sites. Thus, as for the case of perfect state transfer, the 
perfect wave-packet splitting can be a new tool for varied applications. 


I. INTRODUCTION 


The quest for a quantum computer is boosting the develop¬ 
ment and engineering of new sophisticated quantum devices 
that allow us to observe the space-time evolution of its con¬ 
stituents. Indeed, in recent years several experimental groups 
successfully measured the quantum dynamical evolution of 
particles and/or quasi-particle hopping in a lattice IMl. Due 
to the inherent nature of quantum mechanics, the evolution of 
an isolated quantum system is represented by a wavefunction 
t) which describes the probability amplitude of finding a 
particle in position x at time t. Quantum interference can give 
rise to particular structures and patterns in the space-time evo¬ 
lution |0^(x, OP which are known as “quantum carpets” m, 
quantum revivals ifTTIl . or quantum Talbot effect lIT^ . quan¬ 
tum walks (TsHISl, and quantum self-imaging CSl. 

An interesting case is when the wavefunction undergoes a 
revival, namely when after a particular time the shape of the 
initial wave packet is almost perfectly reconstructed. Aside 
from its fundamental implication, revivals occurring into a 
different position, far from the initial one, are particularly im¬ 
portant for connecting and linking distant quantum registers 
(TTlITSl. On the other hand, a lattice of static localized parti¬ 
cles represents an alternative paradigm for quantum commu¬ 
nication where information carriers are not physically moving 
particles but rather collective excitations whose space extent 
is reconstructed at a different position after a certain time. In 
this respect, spin chains represent one of the most viable solu¬ 
tion and there are various protocols to exploit their dynamics 
for transferring states and entanglement between remote sites 
|[T9l l20l . The coherent excitation transfer, or in general the 
wavefunction reconstruction at a certain time, corresponds to 
the phase alignment of the eigenstates entering into the wave 
packet and, as such, can happen only when the energy eigen¬ 
values satisfy certain conditions ||2TJ[22]|. Some models ad¬ 
mitting a perfect |[2T][23l[24]| or almost perfect 125^281 recon¬ 
struction have been explicitly constructed. On the other hand, 
if the phase alignment is only between particular subsets of 
the energy eigenstates, then the wavefunction is split into a 
superposition of copies of the initial wave-packet, each sepa¬ 


rated by a certain distance. This effect is known as fractional 
revival irmirni^. or fractional Talbot effect m. 


In this paper we engineer a chain with nearest neighbor 
interactions to obtain a perfect wave-packet splitting and re¬ 
construction during a ballistic evolution. In other terms, 
if t=U)=f{x) is the shape of the initial wave-packet, at 
the revival time f the wavefunction is if/{x, ^ 

where v is the group velocity defined by the energy eigenval¬ 
ues. While in general the revival time is connected to specific 
algebraic properties of the spectrum and might be very long, 
in our case the splitting happens on a time which is dictated by 
the group velocity of excitations and, as such, scales only lin¬ 
early with the distance. Our method is therefore specifically 
targeted for applications where a smaller operational time is 
particularly important for neglecting the interaction with the 
surrounding environment. Recently, it has been shown that 
the wavefunction of a one-dimensional excitation can be split 
into a transmitted and reflected components by introducing 
localized impurities |[30H32ll . or via suitably designed time- 
dependent control fields 1^ . Here we focus on a different 
strategy aiming at obtaining a perfect fractional revival. 


The generalization of the fractional revival to a many- 
particle setting has many important applications. As far as 
identical particles (bosons/fermions) are concerned, it allows 
one to define a perfect effective beam splitter operation be¬ 
tween distant sites and then to observe multi-particle Han¬ 
bury Brown and Twiss interference effects 1^ . such as per¬ 
fect bunching or anti-bunching. As for spin systems, that the 
perfect fractional revival can be used to generate dynamically 
long-distance entanglement, a topical application which may 
be tested experimentally with current technology ESI ESI. 
Indeed, the use of particle delocalization to generate entan¬ 
glement is particularly evident in a single excitation setup, 
namely when there is a single spin in the state |T) while all 
the other spins are in the state |^). If the wave-packet of 
this excitation is perfectly split and reconstructed in two dis¬ 
tant sites n and m, then the final state of the spins pair {n, m) 
is (ITvDwm + \l'\)nm)l V2, namely a maximally entangled Bell 
state. We show that this reasoning can also be used in a mul¬ 
tiple excitation scenario to dynamically generate a maximal 
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set of Bell pairs in a spin chain setup, and to provide a more 
general version of previous proposals |[J7ll38]| . 

This paper is organized as follows. In section II, we de¬ 
fine the mathematical conditions which allows a particular 
fractional revival, namely the perfect splitting and reconstruc¬ 
tion of an incoming wave packet, and we propose a numeri¬ 
cal algorithm to find suitable Hamiltonians which fulfill these 
conditions. Interesting applications are then analyzed in sec¬ 
tion III in a many-particle setting. In particular, we discuss 
bunching/anti-bunching effects in atoms trapped in an optical 
lattice and the dynamical generation of entanglement in spin 
chains interacting with nearest-neighbor XY couplings. Con¬ 
clusions are drawn in section IV. 


II. PERFECT SPLITTING WITH ENGINEERED 
COUPLINGS 


We study a one-dimensional quantum walk in a lattice 
with nearest-neighbor engineered interactions described by 
the Hamiltonian 


H,= 


ri-i 

Jn\n){n+\ \ + h.c. 

\n=l 


L 


Y^B„\n}{n\, 

n=\ 


( 1 ) 


where \n} represents the state where a particle is in the ^-th 
site, and L is the length of the chain. To find the mathematical 
conditions for a perfect splitting and reconstruction, we first 
focus on the requirements to achieve perfect state transfer. To 
perfectly transfer an excitation from site 1 to site L the co¬ 
efficients Jn and Bn have to satisfy some conditions (see e.g. 
Ref. (221 ). Firstly, the Hamiltonian has to be mirror symmet¬ 
ric, i.e. JL-n=Jn and BL+i-n=Bn for any 1<^<L. The mirror 
symmetry imposes some relations between the eigenvectors of 
the Hamiltonian (241 : if the eigenvalues of Hp are ordered 
such that Ek<Ek+\, then 

Oik = , (2) 


where Hp = OEO^ is the eigenvalue decomposition of Hp. 
The second requirement is that the energy eigenvalues E^ sat¬ 
isfy the relation 


where f is the transmission time and a is an arbitrary phase. 
Here we consider a=0, namely TrH=0. Among the analytic 
solutions of ([^, the simplest one is given by the coupling pat¬ 
tern (231 [391 

yf T = g ^n(L-n), B„ = 0, (4) 

which implements perfect state transfer (PST) at f=LIJ. 
Other solutions can be obtained numerically using inverse 
eigenvalue algorithms EasnisD. If the eigenvalues and 
eigenvectors of Hp satisfy Eqs.(|^ and ([^, then 

k k 

~ ^ On,k^L+l-m,k ~ ^n,L+l-m ? ( 5 ) 

k 


namely an excitation initially located in site m is perfectly 
transferred to site L-m-rl after a time f. 

In a similar fashion, a perfect wave-packet splitting and re¬ 
construction can be obtained when the eigenvalues of Hp sat¬ 
isfy 


e = cos 6 -r /(-l)^ sin 6 , (6) 

for some angle 6. Indeed, by repeating the calculation Q one 
finds {n\e~^^p^* \m)= cos 6 sin 6 Sn,L+i-m ^ namely 

f 

\m) cos^lm) -I- /sin^|L-m-l-l) . (7) 

The eigenvalue relations are one of the main result of this 
paper. By properly choosing 6 it is possible to balance the re¬ 
construction on distant sites, as show in Eq.Q, and for 7r=7r/4 
one obtains the perfect delocalization between distant sites of 
an initially localized wave packet. The coupling pattern to 
satisfy Eq.([^ can be obtained using inverse eigenvalue tech¬ 
niques. Erom the conceptual point of view an inverse eigen¬ 
value problem deals with finding the zeros of the highly non¬ 
linear function f(A)=E(A)-E, where the vector E(A) contains 
the ordered eigenvalues of the Hamiltonian Hp(A) with param¬ 
eters A, and the vector E contains the target spectrum. Among 
the algorithms to find the optimal parameters (42l [43l, the 
most used one relies on the application of the Newton method 
to find the zeros of f(A). The Newton method starts with an 
initial guess and updates it according to the rule (44ll 

/'’)]=/(/”)), (8) 


where the matrix, with elements 






= ,„l Oal">)|t). 

(9) 


is the Jacobian matrix and Hp{A)=0{A)E{A)0{AY is the 
eigenvalue decomposition of Hp{A). The linear system has 
a unique solution provided that is an invertible matrix. This 
in turn implies that the number of parameters have to match 
the number of eigenvalues, i.e. the dimension of the matrix. 

The mirror symmetric Hamiltonian ([T]) has L independent 
parameters, being L the number of sites. Indeed, because of 
the mirror symmetry, when L=2N (being N an integer) there 
are N independent values of Jn and N independent values 
of Bn. On the other hand, when L=2N+\, there are N in¬ 
dependent values of Jn and N-\-\ independent values of Bn. 
We apply inverse eigenvalue techniques to find the coupling 
pattern which allows a perfect balanced splitting of the wave- 
packet. The latter is obtained by imposing the condition 
with 6=7114, so the target eigenvalues are 


LE 

~T 



71 71 

4’4’ 


- - -I- 27T, - -h 271, - 
4 4 


-h 471, . 



( 10 ) 


where, without loss of generality, we have imposed f = LjJ. 
Because f{A) is in general a non-convex function, possibly 
with many local minima, inverse eigenvalue problems are 
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(a) Even length chain 



FIG. 1. Comparison between the perfect wave-packet splitting cou¬ 
plings and the perfect state transfer couplings in Eq.(|^ for 
(a) an even chain and for an (b) odd chain. Only in the latter case the 
perfect splitting requires also the engineering of a field profile 


known to converge only if the initial guess is not too far 
from the ideal set of parameters A for which E(A)=E 1441 . We 
guess that the optimal parameters for a perfect wave packet 
splitting are given by a local perturbation of the fully engi¬ 
neered chain which guarantees perfect state transfer, so we 
use the coupling pattern 0 as an initial condition. 

The algorithm quickly converges to an optimal parameter 
set and hereinafter we called and the obtained opti¬ 
mal couplings and local fields. Surprisingly, we find that for 
even L the algorithm always converges to a solution where 
while for odd L the local fields are different 
from zeros especially near the center of the chain. For exam¬ 
ple, the Hamiltonians for L=5,6 are shown in Appendix A. 
The output of the algorithm is shown in Fig |l(a)| for L=50, 
and in Fig 1 1(b) [ for L=49. As it is clear, both for L even and 
odd the coupling patterns for perfect wave-packet split¬ 
ting are similar to the coupling pattern , in formula Q, 
for perfect state transfer: the only difference being the pres¬ 
ence of few impurities at the center of the chain. Moreover, 
for odd L one requires also the engineering of the local fields 
according to some particular profile. The resulting field 
pattern is constant far from the center of the chain and has a 


particular oscillatory profile near the central sites. 


III. APPLICATIONS 


A. Perfect bunching/anti-bunching in a bosonic lattice 


As a concrete application of the results of the last section 
we consider a model of hopping particle in a one-dimensional 
lattice, described by a Bose-Hubbard Hamiltonian with site 
dependent parameters 


L L L 

Ei — — ^ Jn {^n^n+\ F h.C.^ -h ^ JJnl^n “ 1) “ ^ 

n=l n=l n=l 

( 11 ) 

Here are the tunneling matrix elements, Un is the 

onsite interaction and Bn=B^n^^^ is the chemical potential, Un 
is the boson annihilation operator and nj=a\aj. The Hamil¬ 
tonian O accurately describes cold bosonic atoms in op¬ 
tical lattices US] |46l, and it also models fermions (471 SHI 
and hard-core bosons mg dependently on the onsite interac¬ 
tion values. The tuning of the site dependent coupling con¬ 
stants in O is achieved via addressable optical lattices (501 . 
created projecting an electric field profile via holographic 
masks Elllsll or via micro-mirror devices mil. Initialization 
and read-out of single atoms are achieved exploiting single¬ 
particle addressing techniques dEnsma while magneti¬ 
cally induced Feshbach resonances allow a global control of 
the onsite interaction acting on the collisional coupling con¬ 
stants values (5^ . For instance the non-interacting regime 
has been recently achieved with this technique using Cs 
atoms loaded in a one-dimensional optical lattice (53. 

Thanks to the techniques developed in this paper, the cou¬ 
pling profile produces a splitting of a single particle wave- 
function, which is reconstructed at the transfer time f as 
two copies of the initial wavepacket with probability 1/2 
each. More precisely when the coupling pattern is imple¬ 
mented, the wavefunction of a bosonic atom initially onsite n 
is split by the impurity pattern at the center of the lattice and, 
at the transfer time f, that particle is perfectly delocalized be¬ 
tween two mirror symmetric sites, n and L-n+\. If an another 
particle was in the lattice in position m, after f it would be de¬ 
localized between the sites m and L-m-i-1. When two particles 
are initially in two mirror symmetric sites, i.e. m=L-n-\-\ the 
dynamics generates multi-particle Hanbury Brown and Twiss 
correlations (Ml at f. Indeed, in the free boson case, namely 
when t/„=0, because of the symmetries of the bosonic wave- 
function, after a time f the state becomes 




\^)n\^)m + |0)n|2)m 


\^h), 


( 12 ) 


i.e. the output state consists of a superposition of two bosons 
being in site n and two bosons being in site m=L-n+\. This 
“bunching” effect is the celebrated Hong-Ou-Mandel effect 
(HOM) (23 which has been observed recently exploiting the 
coherent evolution of two particles in a single double-well tun¬ 
neling model (59l . With the results presented in this paper. 


















4 


because of the perfect reconstruction of wave-packets at the 
transfer time, it is possible to achieve a perfect bunching be¬ 
tween arbitrary distant sites of an optical lattice. On the other 
hand, in the strong interacting case, namely in the hard-core 
boson limit f/„=oo, the final state is ))=|1)^|1)„, i.e. there 
is one particle in position n and one particle in position m. 


1. Effect of imperfections in tuning the parameters 


In real systems random noise effect, due to environmental 
variables, and engineering imperfections in the coupling con¬ 
figurations produce deviations from the ideal coupling values 
Eo). The effect of the coupling randomness, even for non¬ 
interacting systems, is to produce a localization of the eigen¬ 
states of the system and consequently to inhibit the state trans¬ 
fer (601. We also mention that recently it has been shown 
El ED that the interaction of bosonic atoms with static 
fermionic impurities, randomly distributed in the lattice, may 
yield a Bose-Hubbard model where the parameters and Bn 
are subject to noise. Given the above, we investigate what de¬ 
gree of imperfections is tolerable in our scheme or, in other 
terms, what is the precision required in tuning the coupling 
strengths according to the desired pattern. 

We firstly include an off-diagonal disorder term (hopping 
disorder) in the Hamiltonian (TT\ as where 

[-6,6] is a uniform random distribution and e is the per¬ 
turbation strength (631 . In Fig. |2(a)| the relative variation 
|AFii |/Fii( 6=0) is shown as function of the degree of disorder 
6. Here Pn=\{^b\^0*)}\^, where is defined in (H) , and 
AF=|Fii(6)-Fii( 6=0)1 represents the deviation of the bunch¬ 
ing probability respect to the ideal case. We also consider the 
effect of diagonal noise Bj=B-\-Jxj with [-rj, rj] in an even 
site chain. The effect of signal noise is shown in figure Fig. 
2(b)| as function of the noise coupling strength 6. As clear 
from Fig. |2(a)| and |2(b)| a power law behavior, under a cer¬ 
tain threshold value of e and rj, characterizes both the devia¬ 
tions due to hopping disorder and due to the diagonal disorder. 
Clearly, an high degree of disorder produces state localization, 
which completely destroys the effect. 


B. Quantum many-particle carpets 

The perfect reconstruction scheme developed in the pre¬ 
vious sections allows generating periodic space-time quan¬ 
tum interference patterns of multi-particles systems known as 
“quantum carpet”. By using the engineered chain with 
and B^n^^^ one in fact expects a regular temporal pattern in 
the evolution: the wave-packets composing the initial state 
are split into two copies, reconstructed into different positions 
after the time f, and then they go back to the initial posi¬ 
tion after a time 2f . On the other hand, during intermediate 
times, quantum interference leads to different behaviors which 
are expected to be susceptible to the particle statistics. To 
show this effect, we study the quantum carpet generated by the 
space-time evolution of the mean occupation number (nfft)}, 



(a)OfF-Diagonal Noise 



(b)Diagonal Noise 

FIG. 2. Relative variation of the bunching probability Pift^f) in 
the non-interacting regime C/„=0, in presence of (a) random hopping 
coupling strength e and (b) random diagonal coupling strength t]. 
Several chain lengths L are considered. 



Site 

FIG. 3. Quantum Carpet: space-time evolution of {n^j(t)) for a L=40 

one-dimensional chain initialized in |0^o)=^Io^3ilO)- We consider 
the free-boson regime C/„=0. The white zones are out of the range 
[0,0.4] which has been chosen for convenience. 


or by the square occupation mean {nj(t)}. The regular inter¬ 
ference pattern of a two particle system is depicted in Fig. 
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FIG. 4. Quantum Carpet due to four particle interference. The initial state \i//o)=alalalal^\0) contains four bosonic particles. We consider the 
space-time evolution of {rijit)) in figure (a,b,c), and the space-time evolution of {n^j{t)) in figures (d,e,f). The chain length is L=14 and several 
value of the onsite interaction Un-U are considered. Here f is the fractional revival time, while 2f is the full revival time. The difference 
between the first and second row is due to bunching/anti-bunching effects. Note the transition from bosonic (t/=0) to fermionic and hard-core 
boson {U-oo) behavior as a function of U. 


where we show the expectation value for two non inter¬ 

acting bosons initially in \i//q)= 01 ^ 0 , 1 ^ 10 } in a one-dimensional 
chain with L=40. 

To highlight more in detail the multi-particle statistical in¬ 
terference effect we consider a system of four particles, ini¬ 
tially in \i//Q)=a2a^^ala[^\0), where L=14. We show in Fig. 
Qa) and|^d) respectively the mean occupation number and 
the quadratic mean occupation number for the non interacting 
case and for the strong interacting case in Fig. |^c) and^f). 
In the boson case bunching effects are observable at f while 
in both cases a perfect reconstruction of the initial wavepacket 
happens at 2f . This is evident more clearly in Fig. where we 
represent the mean occupation number and the quadratic mean 
occupation number of site 2 as function of time. We also take 


into consideration the role of the onsite interaction which af¬ 
fects the perfect reconstruction of a two particle wavepacket. 
It turns out that from the space-time dynamics of {njit)) it 
is not possible to discriminate free evolution (f/=0) from the 
hard-core limit (U=oo), while particle statistics give rise to 
different dynamics for On the other hand, for inter¬ 

mediate values of the onsite interaction the dynamics does 
not lead to a perfect reconstruction of a wavepacket, due to 
scattering effects. This effect is clearly shown in Fig. 0b), 
0e) and in Fig. 0for C/„=l where {r^{t=2f))<{r^{t=Q)) and 
lfl2it=2f))<{n2(J=G)). 
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(a)Mean occupation number 



(b)Quadratic mean occupation number 

FIG. 5. Plot of the (a) mean occupation number and of the (b) 
quadratic mean occupation number of site 2 as function of time for 
several values of the onsite interaction. Here L-IA and the initial 
state is \\l/{0))=a\a\a\a\^\0). Note the transition from bosonic (C/=0) 
to fermionic and hard-core boson (JJ-oo) behavior as a function of 
U. 


C. Perfect generation of entanglement in an XY spin chain 

We now consider a chain of spin-^ magnets described by 
the XY Hamiltonian 

- X “ Tj 

n n 

where a=x,y,z are the Pauli spin operators acting 
on the spin localized in the ^-th site of the chain, and 
(T^=((T^+/cr^)/2. Effective spin-^ systems coupled by the 
Hamiltonian ( p^ with site dependent coupling strengths can 
be obtained in different physical realizations; e.g. in NRM 
using global rotations and suitable field gradients |[64l . with 
atomic ions confined in segmented microtraps l65ll . with neu¬ 
tral atoms trapped into an optical lattice by polarized laser 
beams |[3l|66l, or with superconducting qubits coupled either 
by site dependent capacitors El or inductors |[ 68 ]| . 

The system Hamiltonian ( p3] ) can be mapped to a fermionic 
hopping model via the Jordan-Wigner transformation: the op¬ 
erators Cn= Yl Satisfy canonical anticommutation 

relations and clcm , where H is the hopping 

matrix Q- Every many-body spin state can be obtained by 
applying the annihilation operators to the fully polarized 


state |Q)=|TT • • •)• Therefore, the time evolution of a generic 
initial state can be obtained by expressing the operator in 
the Heisenberg picture 1^ as 

Cn(t) = c„. (14) 

m 

We now show how one can create entanglement between 
two remote mirror symmetric sites by exploiting the perfect 
wave-packet splitting Q. Suppose that, starting from the fully 
polarized state |Q) a particle is flipped in position n\ the ini¬ 
tial state of the system is then c„|Q). When the single-particle 
Hamiltonian implements the transformation 0 , then, thanks 

to Eg. (pi} one has c„|Q) ^ (c„|Q)-f/cL-n+i|0))/V2. There¬ 
fore, going back to the spin picture, after the time f an entan¬ 
gled state between sites n and L-n+\ is generated. 

The above arguments can be generalized in a many-particle 
setting to generate the maximal amount of entangled pairs 
starting from a separable state. Two suitable choices of the 
initial state are 

= (15) 

|^AFM^ = |tiT---n), (16) 

namely the domain-wall state or the anti-ferromagnetic 

state If the system is initialized in either or 

I 0 .AFM) Yqi to evolve under the perfect splitting Hamilto¬ 

nian, then the resulting state after a time f is (ci -f e^^^CL)(c2 + 
e^^^CL-i){c3 -f e^^^CL-2) • • • 1 ^), where at depends on the initial 
state. By carefully dealing with the Jordan-Wigner phase en¬ 
tering into the definition of the operators one can easily find 
that the resulting state corresponds to a state in which every 
pair of qubits sitting in positions n and L-n-\-l is maximally 
entangled. The perfect splitting dynamics thus represents an 
alternative to other methods existing in the literature to gener¬ 
ate nested Bell pairs (371 Ell starting from a separable state. 
However, compared to previous proposals it is more general 
because it allows tuning the number of generated Bell pairs 
by simply choosing the number of flipped spins in the initial 
state. 


IV. CONCLUSIONS 

In this paper we study the wavefunction dynamics of hop¬ 
ping particles and/or quasi-particles in a quantum chain. 
We design the Hamiltonian so that a localized wave packet 
evolves coherently along the chain without dispersion, and at 
particular point is perfectly split into transmitted and reflected 
components which propagate in opposite directions without 
dispersion. When the reflected component reaches the ini¬ 
tial site, its wave packet becomes localized while, at the same 
time, the wave packet of the transmitted component becomes 
localized in a different site of the chain. We devise the exact 
conditions that the Hamiltonian spectrum has to satisfy to al¬ 
low for the perfect splitting and reconstruction. Then we focus 
on some viable Hamiltonians with nearest-neighbor interac¬ 
tions and site-dependent couplings, and we And the coupling 
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pattern which satisfies the perfect splitting condition using in¬ 
verse eigenvalue techniques. 


Besides shedding new light into quantum interference phe¬ 
nomena in one dimension, our results are particularly use¬ 
ful for applications. In this respect, we study atomic lattices 
and obtain perfect Hanbury Brown and Twiss correlations and 
peculiar quantum interference patterns which result in regu¬ 
lar structure in the space-time evolution of the many-particle 
wave function. Moreover, we show that in a spin chain set¬ 
ting, the particle splitting can be used to generate maximally 
entangled states between distant parts. 


We expect that the perfect wavepacket splitting will become 
a general tool for varied applications in controlled quantum 
interference and quantum information processing. 
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Appendix A: Perfect splitting Hamiltonian for L = 5,6 

The Hamiltonian matrices {n\H\m) for perfect balanced 
splitting when L=5 and 6 are respectively (in unit of /): 

f-0.08378 0.6195 0 0 ^ ] 

0.6195 -0.2932 0.6664 0 0 

0 0.6664 0.7540 0.6664 0 

0 0 0.6664 -0.2932 0.6195 

. 0 0 0 0.6195 -0.08378 . 

( 0 0.5999 0 0 0 0 >1 

0.5999 0 0.8279 0 0 0 

0 0.8279 0 0.3927 0 0 

0 0 0.3927 0 0.8279 0 

0 0 0 0.8279 0 0.5999 

. 0 0 0 0 0.5999 0 . 

As shown in Section III, small imperfections parameter tuning 
result in negligible deviations from the ideal dynamics. 
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